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POINT  SOURCE  RADIATION  PATTERN  SYNTHESIS 
BY  ITERATIVE  TECHNIQUES 

1.  Introduction 

Many  different  synthesis  techniques  for  antenna  radiation  patterns 
have  been  developed  over  the  years  [1-6].  Unfortunately,  most  of  these 
techniques  apply  only  to  a restricted  class  of  problems.  As  one  example 
of  the  case  in  point,  usually  these  techniques  either  apply  to  uniformly 
or  to  nonuniformly  spaced  arrays.  The  reason  for  this  distinction  is 
rooted  in  the  fundamental  different  functional  dependence  of  the  array 
radiation  pattern  on  the  elements'  amplitudes  and  spacings.  The  method 
proposed  in  this  paper  is  independent  of  the  functional  dependence  of 
the  array  parameters  and  therefore  treats  them  on  an  equal  basis.  This 
makes  the  method  very  general  and  applicable  to  a truly  large  class  of 
problems.  Furthermore,  the  method  allows  the  introduction  of  a variety 
of  linear  and  nonlinear  constraints  in  the  synthesis,  giving  the  designer 
the  opportunity  to  simplify  the  antenna  construction.  As  an  example, 
the  element  currents  can  be  constrained  to  have  the  same  magnitude  which 
simplifies  the  feeding  system  considerably.  The  method  is  based  on  the 
well-known  steepest  descent  technique  [7],  implemented  in  the  form  of  a 
computer  program.  Many  other  optimization  methods,  which  are  available 
in  the  literature,  could  also  have  been  applied  to  this  problem  with 
various  degrees  of  success  and  complexity  [8-12].  The  steepest  descent 
has  been  chosen  due  to  its  simplicity  in  implementing  constraints  and 
minimal  programming  effort. 
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2.  The  Steepest  Descent  Technique 

In  order  to  explain  this  technique  let  us  assume  that  we  have  a 
function  of  two  variables  z = f(x,y).  This  is  the  equation  of  a sur- 
face in  the  three-dimensional  space  (x,y,z).  Our  problem  is  to  find  the 
position  (xm,ym)  of  the  minimum  z^  on  this  surface.  The  situation  is 
illustrated  in  Fig.  1.  The  way  we  can  proceed  is  to  choose  any  initial 
values  for  (x,y),  say  (x°,y°),  find  the  corresponding  point  z°  on  the 
surface  and,  starting  from  z , "walk"  always  in  the  direction  of  the 
steepest  descent  until  the  bottom  of  the  valley  is  found.  In  mathematical 
language  the  direction  of  the  steepest  descent  is  that  of  the  gradient  of 
f(x,y)  or  Vz.  We  can  then  go  from  (x°,y°)  to  the  next  point  (x^y1)  by 
computing 

x^  = x°  - kV  z°  (1) 

x 

y^  = y°  - kV  z°  (2) 

3 y 

where  k is  the  gain  constant  which  controls  the  convergence  rate  of  the 
method,  V^z°  is  the  x component  of  yz  computed  at  (x°,y°),  etc.  Once 
(x^,y^)  is  found  we  compute  z^  » f(x\y^)  and  compare  it  with  z°  to  see 
if  we  actually  have  proceeded  down  hill.  We  then  continue  in  a similar 
fashion  until  the  bottom  of  the  valley  z is  found.  The  details  of  this 
method  and  more  sophisticated  ones  are  found  in  the  literature  [7-12]. 

Let  us  now  apply  this  technique  to  our  problem.  In  order  to  obtain 
results  which  are  general  enough  let  us  write  our  equations  for  the 
planar  array  case  where  each  element  may  have  a different  radiation  pat- 
tern P^(6),  oriented  arbitrarily  as  shown  in  Fig.  2.  The  pair  (if^.d^ 
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determines  the  position  of  the  ith  element,  the  pair  (I^.a^)  specifies 
the  magnitude  and  phase  of  the  element  excitation  while  y^  determines 
the  orientation  of  the  element  radiation  pattern.  Let  us  see  how  we 
would  particularize  this  general  representation  for  a few  of  the  dif- 
ferent cases  we  have  in  mind.  For  a linear  array  we  could  position 
all  elements  on  the  x axis.  Therefore  = 0°  or  180°  depending  on 
whether  the  element  lies  to  the  right  or  to  the  left  of  the  origin. 

If  all  elements  have  omnidirectional  patterns  y^  may  have  any  arbitrary 
value  (zero  for  example) . If  the  patterns  are  all  pointing  broadside 
(in  the  y direction)  then  y^  = 90°.  In  the  case  of  an  AM  broadcasting 
array,  as  all  towers  have  an  omnidirectional  pattern  P^S)  = constant  = 1, 
y^  = 0 and  ($^,d^)  will  determine  the  position  of  each  element. 

The  radiation  pattern  of  the  planar  array  consisting  of  N excited 
elements  of  Fig.  2 can  be  written  as 

27rd . 

N — J [— T cos  (<p-<p.)  - a ] 

P(<j>)  - l P.  (<t»-y.)  I.e  (3) 

i-1 

Let  Pg(<}>)  be  the  selected  pattern  of  interest.  Usually  only  the 
magnitude  |Pg (<)>)!  of  the  pattern  is  specified.  The  error  between  the 
computed  pattern  |P(<J>)|  and  |Pg(4>)|  for  any  choice  of  the  array  parameters 
is  given  by 

D(<*>)  = | P (<#-)  | -|Pg(<f>)|  (4) 

if  we  choose  M directions  <f>j  , j = l,2 ,3, . . . ,M,  at  which  we  want  to 
minimize  the  error  D(4>) , we  can  define  the  error  over  the  entire 

synthesis  range  as 

M M 

E-  l W(*  )D2($  ) = [ W($  ) [ | P(<t»  )|  - |P  U )|  ]2  (5) 

j-1  J J j-1  J J s 3 

J 


The  reason  for  the  square  in  the  above  equation  is  to  make  sure  that 
when  E is  a minimum  each  component  of  the  sum  is  also  a minimum  and 
therefore  we  have  the  best  approximation  of  | Ps  ( <#> ) | by  |P(<J>)|  in  the 
mean.  W(^)  is  a weighting  function  that  allows  the  synthesis  pre- 
cision to  be  changed  over  certain  ranges  of  $.  Its  use  will  be  ex- 
plained in  Section  3 in  connection  with  specific  examples.  Note  now 
that  E is  a function  of  all  the  array  parameters 


E = E^,  d i,  li,  oi, 


(6) 


where  is  used  to  represent  the  set  of  variables  3 , ...,<j)^ 

and  etc.  Equation  (6)  is  of  the  type  z = f(x,y),the  only  difference 
being  in  the  number  of  variables.  Instead  of  (x,y)  we  have  now  5N 
variables  for  N elements  each  with  5 variables.  E is  therefore  a "surface" 
in  a 5N  dimensional  space.  The  steepest  descent  technique  can  then  be  used 
to  find  the  minimum  E^  on  this  surface.  We  proceed  in  the  usual  way.  We 
start  with  an  initial  guess  for  the  parameters  d?,  I?,  a?,  compute 

the  gradient  at  this  point  VE°  and  then  compute  the  new  values  for  the 
parameters  (jJk  , d^ , 1^ , using  equations  (1)  and  (2)  as 


4 - 4 - W 


4 - 4 - kdV° 
4 - 4 - kiV° 

a3  = a?  - k V E° 
“l—i  a a 


y}  - Y?  - k V E° 
-i  -i  Y Y 


(7) 
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iii  these  equations  V^E  are  the  components  of  VE  corresponding 
to  the  variables  evaluated  at  the  point  (j>?,  d?,  1^,  a?,  x?)  and 
is  called  the  gain  constant  for  the  same  variables.  A similar 
interpretation  applies  to  all  other  equations.  A word  of  caution  should 
be  said  about  the  gain  constants.  First  of  all  we  only  follow  the 
steepest  descent  if  all  k's  of  (7)  are  equal.  Second,  if  the  k's  are 
chosen  to  be  too  small  the  method  will  converge  very  slowly  and  if  the 
k's  are  too  large  the  process  may  even  diverge.  These  two  possibilities 
are  illustrated  in  Figs.  3 and  4 respectively.  As  seen  in  Fig.  3 the 
progress  made  in  each  step  is  very  small  while  in  Fig.  4 the  large  value 
of  k causes  the  successive  steps  to  oscillate  around  the  minimum.  The 
choice  of  k is  very  critical  and  has  been  treated  extensively  in  the 
literature  [7,9,10,11,12].  The  reason  for  the  use  of  different  k's  in 
(7)  is  that  it  allows  the  implementation  of  certain  constraints  very 
easily.  For  example  if  we  desire  to  maintain  the  currents  all  constant 
and  equal  to  the  initial  guess  we  just  have  to  set  k^  = 0. 

Any  types  of  constraints  are  easily  implemented  in  this  method. 

Figure  5 shows  the  case  in  which  the  boundary  CC'  cannot  be  crossed. 

At  every  step  of  the  process  we  check  to  see  if  the  values  of  x and  y 

are  such  that  they  represent  a point  in  the  Forbidden  Region  such  as 

pl . If  this  is  the  case  the  value  of  k is  adjusted  such  that  we  stop 

just  on  the  boundary  at  P^.  If  at  P^  the  gradient  still  points  to  the 

2 

Forbidden  Region  we  are  forced  to  "walk"  on  the  boundary  as  in  P 
otherwise  we  move  to  the  Permissible  Region  and  proceed  in  the  usual  way. 
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Application  to  some  Specific  Synthesis  Problems 


Let  us  now  apply  the  above  ideas  to  some  specific  cases. 

Example  1 - Assume  we  want  to  synthesize  the  pattern  specified 
by  the  eighteen  small  circles  (M  = 18)  in  Fig.  6 for  an  AM  Broadcasting 
Antenna.  Therefore  the  elements  are  omnidirectional  implying  P^S)  = 1 
and  = 0.  We  want  to  use  3 elements  (N  = 3)  and  do  not  have  any 
specific  requirement  in  the  elements  position.  The  synthesis  resulted 
in  an  almost  linear  array  in  the  270°  radial.  The  element  spacing  is 
approximately  A/4  and  the  currents  are  almost  equal.  We  could  have 
performed  the  synthesis  imposing  the  constraint  that  the  currents  be  all 
equal  and  that  the  array  is  linear.  Instead  of  doing  this  let  us  proceed 
to  a more  interesting  example. 

Example  2 - The  desired  radiation  pattern  is  specified  by  the 

twenty-two  small  circles  (M  = 22)  of  Fig.  7.  Let  us  try  to  synthesize 

this  pattern  with  four  elements  (N=4) , each  having  the  radiation  pattern 

2 

shown  in  Fig.  8 (cos  6 with  10%  back  lobe),  and  let  us  assume  that  the 
elements  are  fixed  in  position  and  orientation  as  shown  in  the  insert 
of  Fig.  7.  Therefore  the  variables  in  the  problem  are  only  the  current 
amplitudes  since  we  desire  to  maintain  the  phase  constant  too.  With  so 
many  constraints  we  should  not  expect  to  obtain  a good  fit  as  shown  in 
Fig.  7. 

Example  3 - Let  us  now  use  the  same  data  of  Example  2 but  let  all 
parameters  change  freely.  The  result  is  shown  in  Fig.  9 with  a consequent 
marked  improvement  in  the  fit.  The  array  parameters  as  well  as  the  elements 
are  also  shown  in  Fig.  9. 
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Example  4 - One  of  the  authors  has  developed  the  theory  of  skewed 
arrays  [13,14].  Let  us  try  to  obtain  a pattern  as  smooth  as  possible, 

| P (0) | = 1,  using  four  elements  with  the  radiation  pattern  shown  in 
Fig.  8 but  with  no  back  lobe.  The  elements  are  constrained  to  be  in  the 
radials  0,  90°,  180°,  270°,  the  distance  to  the  origin  being  the  only 
parameter  that  change.  Each  element  has  its  pattern  oriented  in  a 
skew  fashion  as  shown  in  the  insert  of  Fig.  10.  The  amplitudes  are  to 
be  equal  to  one  and  the  phases  equal  to  zero.  The  result  is  shown  by 
the  solid  line  in  Fig.  10  and  represents  an  improvement  over  the  theo- 
retical result  (shown  dashed)  obtained  previously  [14].  The  inter 
element  distance  is  3.64  * Si  5.2X  instead  of  5X  as  theoretically 
predicted . 

Example  5 - The  same  problem  of  Example  4 is  used  here  with  the 
difference  that  now  each  element  has  a pattern  like  Fig.  8,  including 
the  10%  back  lobe.  The  result  is  shown  in  Fig.  11  which  is  still  better 
than  the  theoretical  one  [14],  The  inter  element  separation  is  the  same 
as  that  of  Example  4. 

Example  6 - In  this  example  we  will  try  to  synthesize  the  vertical 
pattern  of  a TV  transmitting  antenna.  The  specified  pattern  is  shown 
in  Fig.  12  by  the  small  circles.  Each  element  has  the  radiation  pattern 
shown  dashed.  The  antenna  feed  system  is  shown  in  the  insert.  The 
following  constraints  are  imposed: 

1 - All  elements  should  have  equal  current 


%*.  ' 
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2 - The  element  spacing  should  be  greater  than  7.5A,  since  each 

element  has  approximately  7A  aperture.  The  two  central  ele- 
ments should  be  8A  apart  to  facilitate  the  feed  system 
construction. 

3 - As  the  elements  are  capacitively  coupled  to  the  line  it  is 

desired  that  the  computed  phases  of  each  element  be  given 
by  its  physical  position  along  the  coaxial  line.  This  is 
therefore  a constraint  relating  the  element  position  to  its 
phase . 

We  would  like  to  take  this  opportunity  to  illustrate  here  how  the 
weighting  function  W can  be  chosen.  For  this  particular  case  W was 
selected  in  the  following  way: 

1 - In  the  main  beam  portion  W = 5 since  we  want  the  main  beam  to 

be  synthesized  with  good  accuracy.  This  is  equivalent  to 
specifying  each  point  in  the  main  beam  five  times. 

2 - In  the  "UP"  region  W was  chosen  equal  to  one  or  zero  depending 

on  whether  the  computed  pattern  exceeds  or  not  the  specified 
one.  This  implies  that  we  do  not  care  about  the  particular 
shape  of  the  pattern  in  this  region  as  long  as  it  is  below  the 
specified  level. 

3 - In  the  "DOWN"  region  we  do  not  want  the  pattern  to  be  below  the 

specified  one  since  this  will  result  in  deficient  station 
coverage.  On  the  other  hand  if  the  synthesized  pattern  is  too 
much  above  the  specified  one  the  gain  will  be  decreased,  there- 
fore we  choose  W to  be  equal  to  two  or  one  half  depending  on 
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whether  the  computed  pattern  is  below  or  above  the  specified 
one.  This  means  that  a deviation  below  the  specification 
weights  four  times  more  than  a deviation  above.  The  result 
is  shown  in  Fig.  12. 

Example  7 - If  in  Example  6 we  relax  the  condition  that  the  currents 
have  to  be  of  the  same  magnitude,  then  a better  fit  is  obtained  as  shown 
in  Fig.  13.  In  this  case  the  power  split  between  the  upper  and  the  lower 
parts  of  the  antenna  is  47:53.  It  would  be  a simple  matter  to  resynthe- 
size this  pattern  with  the  constraint  of  a 50:50  power  split.  All  that 
has  to  be  done  is  to  compute  one  of  the  currents  from  the  other  five 
under  a 50:50  power  split  condition. 

4.  Conclusions  and  Recommendations 

It  has  been  shown  here  that  the  steepest  descent  technique,  when 
implemented  through  a computer  program,  is  a very  powerful  tool  for  the 
synthesis  of  radiation  patterns.  The  method  is  extremely  flexible  and 
can  handle  almost  any  kind  of  linear  and  nonlinear  constraints.  It  pro- 
vides results  that  cannot  be  obtained  by  any  other  synthesis  method  known 
to  the  authors. 

The  mutual  coupling  among  the  antennas  has  been  neglected.  This 
effect  can  be  taken  into  consideration  by  using  an  equivalent  point  source 
representation  of  each  antenna  in  its  actual  environment.  This  will  be 
the  subject  of  a subsequent  report. 

In  order  to  use  this  computer  program,  which  is  described  in 
Appendix  A,  the  user  must  have  some  knowledge  of  the  system  for  which 
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it  was  written  (RADC  H-645)  and  be  versed  in  computer  programming.  The 
reason  for  this  is  that  the  main  purpose  in  the  minds  of  the  authors 
was  to  test  the  technique.  The  next  step,  then,  after  mutual  couplings 
are  taken  into  account,  is  to  make  the  program  much  more  user  oriented. 
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APPENDIX 

Computer  Program  Flow  Chart  and  Typical  Listing 

A flow  chart  of  the  computer  program  used  for  the  solution  of  some 
of  the  examples  presented  in  this  paper  is  shown  in  Fig.  14.  An  expla- 
nation follows  making  reference  to  the  numbered  boxes  of  the  chart: 

1 - Certain  statements  of  the  program  have  to  be  changed  to 
specify  different  constraints.  For  example  we  may  desire  that  the 
amplitudes  be  greater  than  .8.  When  any  amplitude  reaches  a value 
smaller  than  that  it  is  set  equal  to  .8. 

2 - The  weighting  function  W (<p ) is  different  for  each  problem  and 
therefore  must  be  specified  when  a new  problem  is  being  worked  out. 

3 - The  program  always  prints  the  specified  pattern  so  that  a 
check  on  its  accuracy  can  be  made  before  the  synthesis  is  started. 

4 - The  computer  will  ask  the  values  of  the  five  variables  speci- 
fied here. 

If:  JGO  = 1 no  output  is  provided. 

JGO  = 2 the  array  parameters  are  printed. 

JGO  = 3 a plot  of  both  the  specified  and  the  synthesized 
patterns  is  provided. 

JGO  = 4 the  synthesized  pattern  is  printed  in  a table  form 

NUM  is  the  number  of  desired  iterations. 

C1,C2,C3  are  the  gain  factors  for  the  amplitudes,  phases  and 

spacings  respectively. 
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5,6,7  - If  Cl  or  C2  or  C3  is  equal  to  zero  no  iterations  on  the 
corresponding  variable  is  performed.  This  is  another  way  of 
introducing  constraints  in  the  problem.  If  Cl  is  zero,  for 
example,  the  amplitudes  will  be  equal  to  the  initial  guess 
throughout  the  synthesis . 

8 - The  switch  JGO  performs  the  functions  explained  in  4. 

The  program  developed  by  the  authors  consists  of  a series  of 
subroutines  or  functions  that  can  readily  be  changed  to  accomodate  a 
new  problem.  Numerical  computation  of  derivatives  has  been  used 
throughout.  If  one  desires  to  use  the  actual  equations  to  compute  the 
gradient  components , all  that  is  required  is  the  writing  of  a new 
subroutine . 

In  the  end  of  this  Appendix  there  is  a listing  of  a program  which  has 
been  specially  tailored  to  the  case  of  synthesis  of  linear  arrays. 

The  following  notation  was  used  in  the  program: 


M 

N 

NN 

C 

cc 

X(l)  to  X(10) 
X(ll)  to  X(20) 
X(21)  to  X(30) 
TET(I) 

SPEC (I) 

DX(I) 


error  after  each  iteration 

number  of  elements 

number  of  specified  points 

2tt/360 

2tt 

amplitude  of  each  element 
phase  of  each  element 

position  of  each  element  referred  to  the  origin 
specified  angles 
specified  field  strength 

derivatives  of  M with  respect  to  each  element 
in  same  order  as  X(I) . 


The  program  is  composed  of  the  following  parts: 
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Lines  00010  - 00110  Computation  of  Field  Strength 

In  this  subroutine  the  field  strength  is  evaluated.  The  routine 
is  called  for  an  angle  TETX  and  returns  the  value  of  the  real  part  of  the 
field  under  the  name  FRE,  the  imaginary  under  FIM  and  the  absolute  value 
under  RES. 

Lines  00130  - 00270  Computation  of  the  Weighted  Least  Squares 

In  this  part  the  weight  function  W is  defined  and  the  value  of  M 

(called  error  E)  is  computed.  The  f simulation  given  in  the  listing,  is 

that  for  a 4 element  Chebyschev  array  with  20  dB  sidelobe  level. 

Lines  00290  - 00390  Computation  of  the  Gradient 

Numerical  calculation  of  the  partial  derivative  of  M with  respect 

to  each  parameter  which  is  a member  of  the  group  denoted  by  NC. 

Lines  00410  - 00530  Alternative  Computation  of  the  Gradient  for  Constant 

Group 

This  subroutine  is  called,  instead  of  the  above,  when  every  parameter  of 
the  whole  group,  for  example  amplitudes,  is  restricted  to  have  the  same 
value,  although  this  value  has  to  be  iterated. 

Lines  00550  - 00610  Iteration 

New  parameter  setting  according  to  Equation  (7) . Here  restrictions 
on  the  magnitude  of  the  variables  can  be  introduced  as  seen  on  line  00600. 
Lines  00630  - 00690  Alternative  Iteration  for  Constant  Group 

New  Parameter  setting  for  the  special  case  referred  to  above. 

Lines  00710  - 00910  Pattern  Output  Graph 

Use  has  been  made  of  a library  routine  called  PLOT.  Both  patterns, 
the  current  and  the  specified  one,  are  plotted. 


! 
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Lines  01Q90  - 01330  Specification  and  Data  for  the  Main  Program 

This  part  of  the  program  is  concerned  with  format  statements  and 

with  the  input  data. 

Lines  01350  - 01680  Main  Program 

Connects  all  the  subroutines  and  functions  as  suggested  in  the 


flow  chart. 


00010 

00020 

00030 

00040 

00050 

00060 

00070 

00080 

00090 

00100 

00110 

00120 

00130 

00140 

00150 

00160 

00170 

00180 

00190 

00200 

00210 

00220 

00230 

00240 

00250 

00260 

00270 

00280 

00290 

00300 

00310 

00320 

00330 

00340 

00350 

00360 

00370 

00380 

00390 

00400 

00410 

00420 

00430 

00440 

00450 

00460 

00470 

00480 

00490 

00500 

00510 

00520 

00530 


* COMPUTATION  OF  F IELDSTRENGTH 
SUBROUTINE  FUNCT (TETX , FRE ,FIM , RES ) 
COMMON  N,NN,C,CC,X (30) 

FRE=0. ; FIM=C. 

STET=SIN (C*TETX) 

DO  30  1=1, N 

FI=C*X (1+10) +CC*X (1+20) *STET 
FRE=FRE+X ( I ) =COS (FI) 

30  FIM+FIM+X (I) *SIN(FI) 

RES=SQRT (FRE**2+F IM**2) 

RETURN; END 


* COMPUTATION  OF  WEIGHTED  LEAST  SQUARES 
SUBROUTINE  ERR (E) 

COMMON  N,NN,C,CC,X  (30) ,TET(64) ,SPEC(64) 
E=0 


Rl=15. ;R2=-15. 

DO  10  J=1,NN 

CALL  FUNCT (TET (J) , FRE, FIM, RES) 
DIFF=RES-SPEC (J) 

IF  (TET  (J)  .GE.  R1 
IF  (TET(J)  .GE.  R1 
IF  (TET  (J)  .LE.  R2 
IF  (TET (J)  .LE.  R2 
IF  (TET  (J)  .LT.  R1 
10  E=E+W*  (DIFF**2) 

RETURN;  END 


• AND. 

DIFF  .LE.  0.) 

W=0. 

.AND. 

DIFF  .GT.  0.) 

W=l. 

.AND. 

DIFF  .LE.  0.) 

W=0. 

• AND. 

DIFF  .GT.  0.) 

w=l. 

.AND. 

TET(J)  .GT  R2 ) 

W=5 

* COMPUTATION  OF  THE  GRADIENT 
SUBROUTINE  DERIV  (NC,E1) 

COMMON  N,NN,CNCC,X (30) , TET (64) , SPEC (64) ,DX(30) 

CO=l .005 
DO  1 J=1 , N 
K= J +NC 

X (K)=X (K) *CO 
CALL  ERR (E2 ) 

DX (K) = (E2-E1) / (X (K) * (CO-1. ) ) 

1 X (K) =X (K) /CO 
RETURN; END 

* 

* ALTERNATIVE  COMPUTATION  OF  GRADIENT  FOR  CONSTANT  GROUP 
SUBROUTINE  DERIV2  (NC , El ) 

COMMON  N , NN , C , CC , X (30) ,TET(64) ,SPEC (64) ,DX (30) 

X (1)=1.005*X (1) 

DO  1 J=2,N 

1 X (J)=X  (1) 

CALL  ERR (E2 ) 

DX (1) (E2-E1)/ (X (1) * .005 ) 

X (1)=X (1)/1. 005 
DO  2 J=2 , N 
X(J)=X(1) 

2 DX  (J) =DX  (1 ) 

RETURN; END 
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00540  * 

00550  * ITERATION  PROCEDURE 
00560  SUBROUTINE  ITER(NC,CK) 

00570  COMMON  NMNNMCMCCX  (30)  ,TET  (64)  ,SPEC  (64)  ,DX  (3(7) 

00580  DO  1 J=1 , N 

00590  K=J+NC 

00598  X(K)+X(K)-CK*DX(K) 

00600  1 IF  (NC  .EQ.  0.  .AND.X(K)  .LT.  .09)  X(K)=.09 
00610  RETURN; END 
00620  * 

00630  * ALTERNATIVE  ITERATION  FOR  CONSTANT  GROUP 
00640  SUBROUTINE  ITER2(NC,CK) 

00650  COMMON  N,NN  ,C';CC  ,X (30)  , TET  (64)  , SPEC  (64 ) , DX  (30) 
00660  X (1)=X (1) -CK*DX (1) 

00670  DO  1 J=2,N 
00680  1 X (J) =X (1 ) 

00690  RETURN; END 
00700  * 

00710  * PATTERN  OUTPUT  - GRAPH 
00720  SUBROUTINE  GRAPH 

00730  COMMON  N,NN,C,CC,X(30) ,TET(64) ,SPEC(64) 

00740  DIMENSION  PR (2) 

00750  TETX=95. 

00760  NNN=37 

00770  CALL  PLOT (TETX, PR, 1.2, -.2, 2,1, NNN) 

00780  DO  12  J=1 , NNN 
00790  TETX=TETX-5. 

00800  CALL  FUNCT ( TETX , FRE , F IM , RES ) 

00810  PR (2) =0. 

00820  PR (2 ) =RES 

00830  DO  13  1=1, NN 

00840  IF  (TETX-TET (I) ) 13,14,13 

00850  14  PR (2 ) =SPEC (I ) 

00860  CALL  PLOT (TETX , PR ,1.2, -.2, 2,0, NNN ) 

00870  GO  TO  12 
00880  13  CONTINUE 

00890  CALL  PLOT (TETX , PR, 1.2, -.2, 2,0, NNN) 

00900  12  CONTINUE 
00910  RETURN; END 
00920  * 

00930  I PATTERN  OUTPUT  - DATA 
00940  SUBROUTINE  DATA 

00950  COMMON  N , NN ,C ,CC , X (30) , TET  (64 ) , SPEC  (64 ) 

00960  56  FORMAT (I3,6G11.3) 

00970  TETX=92 . 

00980  PRINT: "PATTERN  DATA" 

00990  DO  20  J=l, 50 
01000  TETX=TETX-2 

01010  CALL  FUNCT (TETX , FRE , FIM , RES ) 

01020  IF  (FRE*FIM)  3,4,3 
01030  4 PH=0. 

01040  GO  TO  20 

01050  3 PH=ATAN2 (FIM,RES)/C 

01060  20  PRINT  56,  J, TETX, PH, RES 
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* SPECIFICATION  AND  DATA  FOR  MAIN  PROGRAM 
COMMON  N , NN , C , CC , X (30)  , TET  (64 ) , SPEC  (64)  ,DX  (3?) 

5 0 FORMAT  (8H  ERROR=  ,E11.3) 

51  FORMAT  (8H  ERROR=  , Ell.3,10H  DIFF%=  ,G9.2) 

52  FORMAT (7H  AMPL=  ,6G11.4) 

53  FORMAT (7H  PHAS=  ,6G11.4) 

54  FORMAT  (7H  DAMP=  ,6G11.4) 

55  FORMAT (7H  DPHA=  ,6E1].4) 

56  FORMAT (13, 6G11. 3 ) 

57  FORMAT (-H  ,28H  ENTER:  JOB, ITER# ,C1,C2,C3  ) 

58  FORMAT (7H  DIST=  ,6G11.4) 

59  FORMAT (7H  DDIS=  ,6G11.4) 

N=8; NN=35 

PRINT: "X" 

READ:  (X(I)  ,11=1, N)  , (X  (1+10)  ,1=1, N)  , (X  (1+20)  , 1=1, N) 
PRINT: "TET" 

READ: (TET (I) ,I=1,NN) 

PRINT: "SPEC" 

READ:  (SPEC (I ) ,I=1,NN) 

PRINT: "SPECIFICATIONS" 

DO  5 J=1 , NN 

5  PRINT  56 , J , TET ( J ) , SPEC ( J ) 

PI=3. 14159265 
CC=2 . *PI 
C=CC/360. 


* MAIN  PROGRAM 
DO  2 J=l,30 
2 DX(J)=0. 

CALL  ERR (El) 

PRINT  50, El 

10  PRINT  57 

READ : JGO , NUM , Cl , C2 , C3 

IF  (JGO  .EQ.  99)  GO  TO  99 

IF  (NUM  .EQ.0)  GO  TO  7 

6  DO  1 J=1,NUM 

IF  (Cl  .NE.  0.)  CALL  DERIV(0,E1) 
IF  (C2  .NE.  0.)  CALL  DERIV(10,E1) 


IF  (C3 


0.)  CALL  DERIV (20 , El ) 


IF  (Cl  .NE.  0.)  CALL  ITER  (0,C1) 
IF  (C2  .NE.  0.)  CALL  ITER(10,C2) 
IF  (C3  .NE.  0.)  CALL  ITER(20,C3) 
CALL  ERR (E2 ) 

PRINT  51,  E2 , 100* (E1-E2 ) /E2 

I E1=E2 

7 GO  TO  (10,8,11,18) , JGO 

8 PRINT: "STATUS" 

PRINT  52,  (X  (I) , (=1,N) 

PRINT  53,  (X(I+10) ,I=1,N) 

PRINT  58, (X (1+20) , 1=1, N) 

PRINT  54, (DX (I) , 1=1, N) 

PRINT  55, (DX (1+10) ,1=1, N) 

PRINT  59,  (DX  (1+20) ,1  = 1, N) 

GO  TO  10 

II  CALL  GRAPH 
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01640  GO  TO  10 
01650  18  CALL  DATA 
01660  GO  TO  10 
01670  99  STOP 
01680  END 


I 
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METRIC  SYSTEM 


BASE  UNITS: 

Quantity 


Unit 


Formula 


length 

mass 

time 

electric  current 
thermodynamic  temperature 
amount  of  substance 
luminous  intensity 

SUPPLEMENTARY  UNITS: 

plane  angle 
solid  angle 

DERIVED  UNITS: 
Acceleration 

activity  (of  a radioactive  source) 

angular  acceleration 

angular  velocity 

area 

density 

electric  capacitance 

electrical  conductance 

electric  field  strength 

electric  inductance 

electric  potential  difference 

electric  resistance 

electromotive  force 

energy 

entropy 

force 

frequency 

illuminance 

luminance 

luminous  flux 

magnetic  field  strength 

magnetic  flux 

magnetic  flux  density.  « 
magnetomotive  force 
power 
pressure 

quantity  of  electricity 
quantity  of  heat 
radiant  intensity 
specific  heat 
stress 

thermal  conductivity 
velocity 

viscosity,  dynamic 

viscosity,  kinematic 

voltage 

volume 

wavenumber 

work 


SI  PREFIXES: 


metre 

kilogram 

second 

ampere 

kelvin 

mole 

candela 


radian 

steradian 


metre  per  second  squared 

disintegration  per  second 

radian  per  second  squared 

radian  per  second 

square  metre 

kilogram  per  cubic  metre 

farad 

siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

joule 

joule  per  kelvin 

newton 

hertz 

lux 

candela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

pascal 

coulomb 

joule 

watt  per  steradian 

joule  per  kilogram-kplvin 

pascal 

watt  per  metre-kelvin 
metre  per  second 
pascal-second 
square  metre  per  second 
volt 

cubic  metre 
reciprocal  metre 
joule 


Multiplication  Factors 

i ooo ooo  ooo ooo  = in11 

I 000  000  000  = 10* 

1 000  000  = 10* 

1 000  « 10’ 
100  = 10J 
10  = 10’ 
0.1  = 10-  * 
001  = 10-' 
0 001  » 10-’ 
0 000  001  = 10“ 
0 00(1000  001  » 10- * 
o.ooo  ooo  ooo  ooi  ’ itr" 
oooo  ooo  ooo  (too  ooi  = itr" 
0 000  000  000  000  000  001  *=  10-'* 


kg 


A 

K 


mol 

cd 

rad 

sr 

mis 

(disintegration  )/s 
rad/s 

rads 

m 

kg/m 

F 

A-s/V 

S 

A N 

V/m 

H 

V-s/A 

V 

W/A 

V/A 

V 

W/A 

\ 

N-m 

JIK 

N 

kg-m/s 

Hz 

(cycle  )/s 

lx 

im/m 

cd/m 

Im 

cd-sr 

A/m 

VVb 

V-s 

T 

Wb'm 

A 

* 

^ * — - • 

W 

J/s 

Pa 

N/m 

C 

1 

A-s 

N-m 

W/sr 

|/kg-K 

Pa 

N/m 

W/m-K 

m/s 

Pa-s 

m/s 

V 

W/A 

m 

1 

(wave)/m 

N-m 

Prefix 

SI  Syn 

lera 

T 

«'«<• 

c; 

M 

mega 

kilo 

k 

hecto* 

h 

delta* 

da 

dec!" 

d 

centl* 

c 

mllll 

m 

micro 

nano 

pico 

femlo 

n 

r 

alto 

a 

